Various correlations in a Heisenberg XX Z spin chain both in thermal equilibrium and under the 

intrinsic decoherence 
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In this paper we discuss various correlations measured by the concurrence (C), classical correlation (CC), 
quantum discord (QD), and geometric measure of discord (GMD) in a two-qubit Heisenberg XX Z spin chain 
in the presence of external magnetic field and Dzyaloshinskii-Moriya (DM) anisotropic antisymmetric interac- 
tion. Based on the analytically derived expressions for the correlations for the cases of thermal equilibrium and 
the inclusion of intrinsic decoherence, we discuss and compare the effects of various system parameters on the 
correlations in different cases. The results show that the anisotropy ,l z is considerably crucial for the corre- 
lations in thermal equilibrium at zero temperature limit but ineffective under the consideration of the intrinsic 
decoherence, and these quantities decrease as temperature T rises on the whole. Besides, J turned out to be 
constructive, but B be detrimental in the manipulation and control of various quantities both in thermal equilib- 
rium and under the intrinsic decoherence which can be avoided by tuning other system parameters, while D is 
constructive in thermal equilibrium, but destructive in the case of intrinsic decoherence in general. In addition, 
for the initial state |*I/i(0)) = -^=(|01) + |10)), all the correlations except the CC, exhibit a damping oscillation 

to a stable value larger than zero following the time, while for the initial state |$2(0)) = -^=(j00) + |11)), 
all the correlations monotonously decrease, but CC still remains maximum. Moreover, there is not a definite 
ordering of these quantities in thermal equilibrium, whereas there is a descending order of the CC, C, GMD and 
QD under the intrinsic decoherence with a nonnull B when the initial state is | \l/2 (0) } . 

PACS numbers: 03.67.-a, 03.65.Ud, 75.10.Pq 



I. INTRODUCTION 

The most fascinating nonlocal correlation feature of quan- 
tum mechanics is the quantum entanglement, generally con- 
sidered as an essential resource for the quantum information 
processin gjlll that provides the possibility of quantum tele- 
portation [20, quantum dense coding |3], and quantum cryp- 
tographic key distribution [4], etc. Ever since the foundation 
of the quantum information science, the quantification of the 
entangled states has been one of the most fundamental and 
substantial tasks. As is well known, quantum states have 
been subdivided into entangled states and separable (nonen- 
tangled) states. However, recent research reveals that entan- 
glement doesn't provide all aspects of quantum correlations 
which arise from the noncommutativity of operators repre- 
senting states, observables, and measurements JH]. Quantum 
states display other nonlocal correlations not present in the 
classical counterpart, such as the so-called quantum discord 
(QD) that is intimately relevant to local measurement which 
accounts for all nonclassical correlations originally introduced 
by Ollivier and Zurek @]. QD is defined by the distinction 
between the two quantum extensions of the classical mutual 
information defined equivalently in two classical ways and 
shown to be nonzero both theoretically JH 0] and experimen- 
tally |01 for some separable states which may be utilized to 
speed up some tasks over their classical counterparts. In ad- 
dition, QD is responsible for the quantum computational ef- 
ficiency of deterministic quantum computation with one pure 
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qubit in Ref. 10 Si albeit in the absence of entanglement. 
Therefore, it is imperative and desirable to study QD aim- 
ing at understanding well the relationship between QD and 
other correlation indicators and also that among the total cor- 
relations, genuinely classical correlations and purely quantum 
ones. 

Recently, the QD has been intensively investigated in the 
literature both theoretically Ir9t- t32ll and experimentally Ir7[ l33ll . 
Generally, it is somewhat difficult to calculate QD and the 
analytical solutions can hardly be obtained exce pt f or some 
particular cases, such as the so-called X states IU2I1 . Some 
researches show that QD, concurrence (C) and classical cor- 
relation (CC) are respectively independent measures of corre- 
lations with no simple relative ordering and QD is more prac- 
tical than entanglement 10. Quite rencently, Dakic et al l27ll 
have introduced an easily analytically computable quantity, 
geometric measure of discord (GMD), and given a necessary 
and sufficient condition for the existence of nonzero QD for 
any dimensional bipartite states. Moreover, the dynamical be- 
havior of QD in terms of decoherence J 2*1 32^ 34j - [36ll in both 
Markovian [11] and Non-Markovian cases is also 

taken into account. 

In previous studies, the influence of intrinsic (phase) de- 
coherence, a virtually unavoidable effect caused by the inter- 
action of the system with the surrounding environment, on 
the dynamics of various correlations (C, QD, GMD and CC) 
using a Heisenberg spin chain as a quantum channel has not 
been considered. Also the effect of Dzyaloshinskii-Moriya 
(DM) interaction, which is introduced via the extension of the 
Anderson superexchange interaction theory by including the 
spin-orbit coupling effect, on these correlations and the com- 
parison of these quantities have been rarely reported in the 
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literature. On the other hand, due to the good integrability and 
scalability, the solid state systems 139I - I46I1 have gained great 
attention. Particularly, the Heisenberg spin chains, as the nat- 
ural candidates for the realization of the entanglement showed 
some substantial advantages compared with the other physical 
systems I47l - l5lll . In addition, by suitable coding, the Heisen- 
berg interaction alone can be used for quantum computation 
Il52l - D4ll . To this end, in this paper we investigate in detail both 
the thermal equilibrium and dynamical behaviors of various 
correlations in a Heisenberg XX Z spin-A chain. We find that 
the anisotropy J z is considerably crucial for these quantities 
in thermal equilibrium at zero temperature limit but ineffec- 
tive under the consideration of the intrinsic decoherence, and 
these quantities decrease as temperature T rises on the whole. 
Besides, J and D contribute equivalently to these quantities 
and turn out to be the most efficient controlling parameters. J 
plays a constructive role and B a detrimental role in the ma- 
nipulation and control of these quantities both in thermal equi- 
librium and under the consideration of intrinsic decoherence, 
while D is constructive in thermal equilibrium, but becomes to 
be destructive in the decoherent time evolution process in gen- 
eral. Albeit B is detrimental for these quantities, it remains 
worthy of being studied since the inclusion of it is in practical 
need such as the nuclear magnetic resonance quantum com- 
puting and the superconducting quantum computing. Further- 
more, when the initial state is |*i(0)) = 775 (|01) + |10)), all 
the correlations except the CC, exhibit a damping oscillation 
to a stable value larger than zero following the time, while for 
the initial state |* 2 (0)) = ^(|00) + 1 11)), all the correla- 
tions monotonously decrease, but CC still remains maximum. 
Moreover, there is not a definite ordering of these quantities 
in thermal equilibrium, whereas there is a descending order of 
the CC, C, GMD and QD under the intrinsic decoherence with 
a non-zero B when the initial state is |^2(0)}. 

This paper is organized as follows. In section II, we present 
a brief overview of various correlation measured quantities. 
Next we study the thermal correlations in a two-qubit Heisen- 
berg XX Z spin chain with DM interaction in the presence of 
an external magnetic field along the z-axis in section III. Sub- 
sequently, we turn to the influence of intrinsic decoherence on 
various quantities in section IV. Finally we conclude the paper 
in section V. 



IL CORRELATION MEASURES FOR BIPARTITE 
SYSTEM 

Firstly we give a brief overview of various correlation mea- 
sures. Given a bipartite quantum state pab in a composite 
Hilbert space % — Ha <S> T~Lb, the concurrence |55]| as an 
indicator for entanglement between the two qubits is 

C(pab) = max{Ai - A 2 - A 3 - A 4 ,0}, (1) 

where Aj(i = 1, 2, 3, 4) are the square roots of the eigenvalues 
of the "spin-flipped" density operator R = pp = p(a\ <g> 
a V)P*W\ ® ^V) m descending order. a v is the Pauli matrix 
and p* denotes the complex conjugation of the matrix p in the 



standard basis {|00), |01), 1 10) , 1 11)}. For the density matrix 
in the X form, an alternative equivalent expression is given by 

C(pab) =max{C 1 ,C 2 ,0}, (2) 

where Gy = 2(|p 4 i| - VP33P22) and C 2 = 2(|p 32 | - 

In classical information theory 115611 . the total amount of cor- 
relations between two systems A and B can be represented by 
the classical mutual information I(A; B) = H(A) + H(B) — 
H(A, B), where H(A) (H(B), H(A, B)) is the Shannon en- 
tropy H = — pi log 2 Pi , Pi representing the probability of 
an event i associated with A (B, AB). By virtue of the Bayes 
rule, the mutual information can be rewritten as J (A; B) = 
H{A) - H(A\B), in which H(A\B) = H(A, B) - H{B) 
is the classical conditional entropy employed to quantify the 
ignorance of the state of A when one knows the state of B. 
Despite the equivalency of the two expressions in the classi- 
cal case, the quantum versions of the two are not equivalent 
anymore. In the generalized quantum version, the classical 
probability distributions are replaced by the density operator 
p and the Shannon entropy by the von Neumann entropy ll57ll 
S{p)=— Tr(plog 2 p). Accordingly, the quantum version of the 
two mutual information expressions can be obtained as 

1(pa;pb) = S(p A ) + S(p B ) - S(pab), (3) 

where pa(b) = Tib(A) (pab) is the reduced density matrix of 
the subsystem A(B) by tracing out the subsystem B(A). The 
quantum generalization of the conditional entropy is not the 
simply replacement of Shannon entropy with von Neumann 
entropy, but through the process of projective measurement 
on the subsystem B by a set of complete projectors B k , with 
the outcomes labeled by k, then the conditional density matrix 
p k becomes 

Pk = — (lA®B k )p(l A ®B k ), (4) 
Pk 

which is the locally post-measurement state of the subsystem 
B after obtaining the outcome k on the subsystem A with the 
probability 

Pk =Tr[(l A ®B k )p{l A ®B k )}, (5) 

where I4 is the identity operator on the subsystem A. The 
proj ectors B k can be parameterized as B k = V\k)(k\V\k — 
0, 1 and the transform matrix V € U{2) JT^l is 

( cos 9 sin 9 \ 

\ e^sine -cos6 ) ' ( ' 

Then the conditional von Neumann entropy (quantum condi- 
tional entropy) and quantum extension of the mutual informa- 
tion can be defined as Jfl 

S(p\{B k })=^2 Pk S( Pk ), (7) 
fe 

J( PAB \{B k }) = S(p A ) - S(p\{B k }). (8) 
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Following the definition of the CC in Ref.dH] 

CC(p AB ) = sup J(p AB \{B k }), (9) 

{B k } 

then QD defined by the difference between the quantum mu- 
tual information I(pab) and the CC(pab) is given by 

QV{pab) = Arab) - CC(p AB ). (10) 

If we denote S min (p AB ) = mi MJ3k} S(p AB \{B k }), then a 
variant expression of QD reads IU5I1 

QD(PAB) = S(pb) - S(p AB ) + SnaaipAB)- (H) 

It is usually difficult to get the analytical expression of QD 
except for some special cases, thus another correlation mea- 
sure, GMD, is introduced by Dakic et at to simplify the 
computation. It is defined as 

V G {p AB ) = mm\\p-x\\ 2 , (12) 

x 

where the minimum is over the set of zero-discord states and 
the geometric quantity 

\\P-X\\ 2 :=MP-X) 2 

is the square of Hilbert-Schmidt norm of Hermitian operators. 
For any two-qubit state in the so-called Bloch basis 

1 3 1 / 3 

p = - ^ Rij^i ® Vj = J ( lA ® IB + ^(XiCTj <g> 1 B 
i,j=Q ^ i=l 

+ yil B <8> (Ti) + ^ tijCJi (8 (Tj J , (13) 

i,j=l ' 

where = Tr[p(oi <g> a 3 )], <r = hy.2, 0$ = 1,2,3) 
are the Pauli matrices, x = {xi}, y = {j/j} are the three- 
dimensional Bloch vectors associated with subsystems A, B, 
and tij denote the elements in the correlation matrix T. Then, 
a variant expression of GMD is given by 

V G { PAB ) = \{\\yf\\ + \\T\\ 2 - fcmax), (14) 

where fc max is the largest eigenvalue of the matrix yy^ + T T T 
(in the case of measurement on the subsystem A, one needs 
to replace y with x and T T T with TT T ). An alternative for- 
mulation for GMD has been provided in Il6tl . Note that its 
maximum value is | for two-qubit states, so it is appropriate 
to consider 2Dq as a measure of GMD hereafter in order to 
compare with other correlation measures i23tl . 



III. THE CORRELATIONS IN A HEISENBERG XX Z SPIN 
CHAIN IN THERMAL EQUILIBRIUM 

The physical system we discuss here is a two-qubit Heisen- 
berg XX Z spin chain with the DM anisotropic antisymmetric 



interaction l58l l59ll under the external magnetic field and its 
Hamiltonian is written as 

H = \[J (of of + of at) + J z a\al +B(a z 1 + a 2 2 ) 

+ I?(ff?o?-<rM)]. (15) 

where J and J z are the real coupling constants, a\ (i = 
1,2; j = x, y, z) are the Pauli matrices. B and D are respec- 
tively the z-component of the external magnetic field and DM 
interaction. We are working in units, so that all parameters 
are dimensionless. The state of a typical solid-state system at 
thermal equilibrium in temperature T (canonical ensemble) is 
p{T) = e~kr JZ 1 with Z = Tr(e~"Er) the partition function 
and k the Boltzmann constant. Usually we work with natural 
unit system H = k = 1 for simplicity and henceforth. 

In the first instance, we derive the analytical expressions for 
various correlation measured quantities. After some straight- 
forward algebra, one can readily obtain the thermal concur- 
rence 

C(p AB ) =max{i(2e*sinh^- e -*),0}, (16) 

I = 7t _ 2B + .J Z , 2B 

where p = vJ 2 + D 2 and Z = e ar + e t + 

2e ~t ' cosh Next, in order to gain the thermal CC and 
thermal QD, one needs to calculate the quantum conditional 
entropy and minimize it over all possible projective measure- 
ments which is the most difficult part. As per Eqs. ©-(0, 
after the minimization of the quantum conditional entropy 
(i.e., to set the derivative of the quantum conditional entropy 
S(p\{B k }) with respect to angels 9 and <fi to be zero), one can 
find that the quantum conditional entropy is independent of 
angle 0, and reaches its minimum value when 9 = t 2 "'^ 1 )^ 
(m € Z) for — J < J z < J in the absence of B and D. It 
is independent of 9 when J z = ± J (i.e., the XXX model, 
this can be explained from the physical respective as the sys- 
tem is isotropic). Otherwise its minimum value is reached as 
9 = (m e Z). As B and D are introduced, the range 
of J z {—J < Jz < J), in which the quantum conditional en- 
tropy reaches its minimum value when 9 = (2m + 1)7r (m e Z), 
are broadened. Moreover, the effect of B on the dependence 
of the minimum quantum conditional entropy with respect 
to angel 9 is different for the antiferromagnetic (AFM) case 
( J > 0) and ferromagetic (FM) case ( J < 0) with D fixed. 
The effect for the AFM case is more significant compared with 
the case of FM, in which the range of J z {—J < J z < J) 
widens slightly. While the effect of D on the dependence 
mentioned above are the same for both the AFM and FM cases 
when B is fixed. In addition, the range reduces with the rise 
of the temperature. Thereby, the CC and QD are obtained re- 
spectively as 

CC(pab) = S(p A ) - min{Ai, A 2 }, (17) 



QV{p AB ) = S(pb) - S( PAB ) + min{Ai, A 2 }, (18) 
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where 



S(pa) = S(p B ) = -A- log 2 A_ - A+ log 2 A+, 



1 J Z ±2B U 

A± = -(e 2r +e 2r C osri-), 



1 J Z ±2B 1 J^±2 



r/1,2 = 7?e , 773,4 = -e ^ 



Ai = — [In 4 - (1 - 5) ln(l - 5) - (1 + 5) ln(l + 5)], 
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and 



2V(e^-l) 2 + 4e H ^sinh 2 f- 

~ 2B ~ B + J z ' ~ : 

1 + e t + 2e t cosh ^ 



A 2 = ( Vl + w)(£_ + £+) + (m + w)(C- + C+), 



1 ± i/i 1 ± 1 ± i/ 2 , 1 ± u 2 . 



1 i u ^(e^coshf -l) 2 

LJ = —ear cosh ; p 12 



Z 



T 



1 + e t cosh 



Finally, according to Eqs. ([T3T i and (fl4l i. thermal GMD can be 
written as 



2V G ( PAB ) = - - max{ri, T 2 }, 



(19) 



in which 

2e~^cosh^- 



4 cosh ^ cosh | + e ^ (3 cosh 



1) 



4e" 



[(cosh y 



Z 2 

coshf ) 2 +sinh 2 f ] 



Z 2 



4e^ sinh ^ 

r 2 - - 2 . 

In the second instance, we concentrate on the numerical 
analysis of the dependence of various correlations on the dif- 
ferent tunable system parameters at length. Figure Q]plots the 
behavior of various quantities versus temperature T for differ- 
ent isotropy J in the absence of external magnetic field and 
DM interaction (B = D = 0). These quantities are invari- 
ant under the substitutions J — >• — J and D —> — D as well 
as being contributed equally by J and D in thermal equilib- 
rium since J and D only appear in the term p = \/ J 2 + D 2 . 
Without loss of generality, we restrict our attention to the case 
of J > 0. The figure clearly shows that, when J < \J Z \, 
QD begins at zero and increases to a certain value as T rises, 
then decreases with the further rise of T until reaching the 
critical temperature, at which QD vanishes. In addition, this 
phenomenon can only occur for the case of an appropriate J z 



FIG. 1 : (Color online) Various correlations are plotted as a function 
of temperature T for different isotropy J with J z — —0.5 in the 
absence of B and D, in which J = 0.3, 0.5, 1 corresponds to dotted, 
dashed and solid line respectively. QD, CC, C and GMD corresponds 
respecitvely to the red, orange, blue and black line. The inset is the 
uncombined plot. 




FIG. 2: (Color online) Various correlations versus the anisotropy J z 
in the absence of B and D, with J = 1, T = 0.5. The red, orange, 
blue and black line indicates QD, CC, C and GMD respectively. 



in the negative region (i.e. J z < 0). This is also true for GMD 
but not for C and CC. Note that C is always zero in this case, 
while CC starts at maximum and decreases with T. When 
J \Jz\, QD starts at a definite value (for J = \J Z \) and at 
the maximum value (for J > \J Z \), then decreases with T, 
which is also valid for GMD and C except that C is still zero 
when J — \J Z \. However, CC decreases to a certain value 
when J = \J Z \ and increases immediately to maximum for 
J > \J Z \ for low temperatures in the vicinity of zero. The 
results reveal that the quantum phase transition (QPT) occurs 
at J = \J Z \. We should also note that QD is always larger 
than GMD in this case. But there is no definite ordering of 
these measures that are dependent on various system parame- 
ters. Moreover, the critical temperature can be elevated by the 
larger absolute value of the isotropy parameter J. 

Subsequently, in order to demonstrate the effect of J z on 
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FIG. 3: (Color online) Various correlations are plotted as a function 
of temperature T for different magnetic field B with J z — 1, D = 0, 
where B — 0,2,4 corresponds respectively to the solid, dashed and 
dotted line. The red, orange, blue and black line indicates QD, CC, 
C and GMD respectively. The inset is the uncombined plot. 




FIG. 4: (Color online) The QD are plotted as a function of B and 
D (or B and J, or B and J z ) at temperature T = 0.5 with other 
parameters fixed at the value 0. 1 . The upper surface plot represents 
two completely overlapped surfaces QD(B, D) and QD(B, J), and 
the lower one corresponds to QD(B, J 2 ). 



various quantities, we plot Fig. [2] from which one can see that 
both QD and GMD are zero when J z is small in the negative 
region, and start to increase with the increasing of J z at the 
critical point until J z = —J= —1, at which they undergo the 
so-called sudden change. With the further increase of J z , they 
decrease slightly and then continue to increase until reaches 
two different stable values (near the maximum). Note that 
the sudden change also occurs at J z = J = 1 for both of 
them. However the above process is not true for C and CC. C 
is always null in the region J z < — 1 and begins to increase 
abruptly at J z = —1 until to a stable value as J z increases 
and the sudden change does not happen for C. As for CC, it 
decreases from the maximum value to the minimum value as 
J z increases until J z = — 1 and undergoes the sudden change 
twice at J z = ± J, then finally revives to the maximum value 
with the further increase of J z . Thereby, by comparing with 
the result in Ref. Il3ll that QD can signal a QPT at finite tem- 
perature while C can't, we can conclude that not only QD (or 
GMD), but also CC can detect the critical points of QPT. 

Figures [3] and [5] are plotted to exhibit the ffects of B and D 
on various quantities as the external magnetic field and DM 
interaction are introduced. One can clearly see that B plays a 
destructive role in the manipulation and control of these corre- 
lations from Fig. [3] But it still deserves to be studied since the 
introduction of it is sometime in practical need such as the nu- 
clear magnetic resonance quantum computing and the super- 
conducting quantum computing. In particular, the destructive 
effect of B can be compensated through adjusting other tun- 
able system parameters, say J and D for instance. In order to 
compare the efficiency of the parameters against the detrimen- 
tal B, we plot the QD as a function of B and D (or B and J, or 
B and J z ) in Fig. [4] The figure shows that the effects of J and 
D on QD are exactly equivalent and can tune QD to the maxi- 
mal value as long as they are large enough, whereas the effect 
of J z is comparatively weaker (see the lower surface plot in 
Fig. |4|i in compensating the detrimental influence of B on QD. 



The behavior of other three quantities versus these system pa- 
rameters are similar (plots omitted). Furthermore, we should 
note that the difference of these quantities becomes larger as 
B increases until the critical value, at which these quantities 
start at zero and increase as temperature rises (i.e., the value 
of any quantity at zero temperature limit transits from nonnull 
to null at the critical magnetic field or vice versa). With the 
further increase of B, such difference vanishes. Moreover, as 
mentioned above in Fig. [T] QD increases with T when J z 
is in the negative region. Such negative-only region can be 
widened to the positive one by the inclusion of a strong B, 
and also such characteristic is awarded to CC and C. Besides, 
Fig. [5] shows that DM interaction is constructive for various 
quantities. We should note that the difference of QD and CC 
becomes smaller as D increases and the two curves overlap 
each other when D is large enough, while the difference of 
QD and C becomes larger with D. However, QD and GMD 
almost overlap each other all the time and the influence of D 
on their negligible distinction is very slight. In addition, D 
as well as J turn out to be the most efficient parameters in 
increasing various correlations as well as the critical tempera- 
ture. 



IV. THE CORRELATIONS UNDER THE INTRINSIC 
DECOHERENCE 

Now we consider the influence of intrinsic decoherence, 
proposed by Milburn |6()| with the assumption that a system 
does not evolve continuously under unitary transformation for 
sufficiently short time steps, on various correlations. The mas- 
ter equation describing the intrinsic decoherence can be for- 
mulated as 

^- = -i{H,p)-±[H,[H,p(t)}], (20) 
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FIG. 5: (Color online) Various correlations are plotted as a function 
of temperature T for different D with J z = 1, B = 0, where D = 
0, 2, 4 corresponds respectively to the solid, dashed and dotted line. 
The red, orange, blue and black line indicates QD, CC, C and GMD 
respectively. The inset is the uncombined plot. 



where 7 is the phase (intrinsic) decoherence rate. The formal 
solution of the above equation is given by i6lll 



lk 

p{t) = Y,iA Mk{t)mM]k{t) > 

k=0 



(21) 



where p(0) is the density operator of the initial system and 
M k (t) is defined by 



(22) 



M k {t) = H k cxp(-Lfft) cxp ( - —H 



Then the time evolution of the density operator for the Heisen- 
berg XX Z spin system mentioned in the above section can be 
expressed by 



P(*) 



^exp 

mn 



it 



(E m — E n ) 2 — i(E m — E n )t 



X (lp m \p{0)\lpn)\4'm)(lpn 



(23) 



where E m ^ n and V'm.n are the eigenvalues and the correspond- 
ing eigenvectors of the Hamiltonian respectively. 

Firstly, we assume that the initial state of the sys- 
tem is |*i(0)) = ^(|01) + 1 10>) and then pi(0) = 
|^ r i(0))(*i(0)|. Secondly, we consider another initial state 
|*2(0)) = 75(100) + |11)). AspertheEqs. © - Q3) and 
( 1231 , after some algebras, one can obtain the analytical ex- 
pressions for the dynamical behavior of various correlations 
as 



d(t) = trWJ 2 + L> 2 e- 4 ^' cos(2/ii), 
2 4 

QD 1 (t) = - OLi l0g 2 0>i +Y. a i l °&2 a ii 
i—1 i—3 
2 



(24) 



cc 1 (t) = -E"iio g5 



2D G1 {t) = [2J 2 + D 2 e~ 4 ^ t*(1 + cos 2 (4^))] , 
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FIG. 6: (Color online) Time evolution of various correlations are 
plotted as a function of J (with D = 0A),D (with J = 0) and B 
respectively at time t = 1 with fixed decoherence rate 7 = 1, in 
which the orange, blue, black and red line corresponds to CC, C, 
GMD and QD respectively. 



for the initial state |*i(0)) and 
C 2 (f) 



QX> 2 (*) = l+EAlog 3 ft 
<7C 2 (t) = l, 



(25) 



the initial state 



,* 2 (0)), where 
sin(2/i,t)l, 
J 2 e 2 ^^), and 



«1,2 



1 ± 



ifl±e- 2B V) 



for 

M -i e -2 M 2 -ri 1 /D2"- l . j-e-^ ( .p lupi 2 
From the analytical expressions, one can readily see that the 
time evolution of these quantities are independent of J z in 
both cases under consideration. Also, they are independent 
of B in the case that the initial state is |*i(0)), whereas 
independent of J and D in the case that |*2(0)) is chosen as 
the initial state. 

In what follows, we are dedicated to the numerical results. 
In Fig. [6] we plot the time evolution of various quantities 
versus J, D and B respectively with other parameters fixed. 
From the figure one can see that the effects of J and B on 
the time evolution behavior of various quantities are similar 
to that in thermal equilibrium except that CC is always max- 
imal and independent of B. However, the effect of D on the 
dynamics of various correlations are notably different from its 
effect in thermal equilibrium. Note that D is completely detri- 
mental for QD and GMD and nearly destructive for C except 
for a process of sudden death and slight revival, while almost 
beneficial for CC apart from a so-called regrowth process that 
CC decreases to the minimum value (not zero) as D increases 
and then increases to the maximal value. 

Finally, Fig. [7] is plotted for the two different initial states 
in order to observe the effects of pure phase decoherence rate 
7 on the dynamics of various quantities. Before giving the nu- 
merical analysis, we should clarify that the orange-solid line 
and red-solid line in (a) as well as the blue-dotted line and 
black-dashed line in (b) are fully overlapped. Figure [7J a ) de- 
picts that the dynamics of these quantities oscillate with time 
t periodically with the same periodicity and the amplitudes 
of C, QD and GMD decay gradually to a stable value after a 
long time evolution as intrinsic decoherence is taken into ac- 
count, while that of CC, on the contrary, is enhanced with the 
increase of 7. The larger 7 leads to the faster decay (or pro- 



7 



1.00 




FIG. 7: (Color online) Time evolution of various correlations are 
plotted as a function of time t for different decoherence rate 7 with 
two different initial states, where 7 = 0,0.1,0.2,0.3 corresponds 
respectively to the solid, dashed, dotted and dot-dashed line in (a), 
7 = 0.1, 1,2 corresponds to the solid, dashed and dotted line respec- 
tively in (b). The orange, blue, black and red line represents CC, C, 
GMD and QD respectively. 



motion for CC) in a short time. Furthermore, when J = and 
7 is very small, these quantities undergo the sudden death and 
revival periodically. Figure |7Jb) shows that CC is maximal 
all the time but other quantities dissipate monotonously and 
disappear eventually as time goes to infinity so that we can 
conclude that CC is robust against these tunable parameters 
whereas QD is most sensitive. Moreover, there is a descend- 
ing order of CC, C, GMD and QD when the initial state is 
I \E r 2 (0)) with a nonnull B (when B = 0, these quantities are 
all maximal). 



mal equilibrium and under the intrinsic decoherence cases. 
We have obtained analytical expressions for these correlations 
for both cases and discussed their behaviors following vari- 
ous system parameters at length. The results show that the 
isotropy parameter J plays a constructive role in the manip- 
ulation and control of various correlations and the anisotropy 
J z is considerably crucial for these quantities in thermal equi- 
librium at zero temperature limit but ineffective under the con- 
sideration of the intrinsic decoherence. When J z < J < —J z 
(J z is negative) in the absence of B and D, QD and GMD 
start at zero and increase as T rises to a certain value, then 
decrease, while C is zero and CC declines from the maximal 
value with the rise of T. For J = | J z \, all the quantities start 
at a certain value and decrease with T except for C which is 
still zero. When \ J\ > — J 2 , all quantities decrease starting 
from maximum with T, Therefore, J z = ± J are the QPT 
points which are signaled not only by QD and GMD, but also 
by CC. As B and D are introduced to the system, the range 
of J z , in which the quantities start at zero and increase with 
T to a certain value then decrease, are widened to the posi- 
tive region. The inclusion of B turns out to be destructive, 
nevertheless it still deserves to be studied for its practical ap- 
plication in some implementations such as the nuclear mag- 
netic resonance quantum computing and the superconducting 
quantum computing. D plays a constructive role and the ef- 
fects of J and D on the correlations are exactly equivalent 
and they turn out to be the most efficient in compensating the 
detrimental influence of B. Furthermore, the difference of 
various quantities becomes larger with the enhancement of B 
until the critical point, after which it minifies. When the in- 
trinsic decoherence is taken into account, the effect of J and 
B are similar with that in thermal equilibrium, but D becomes 
to be destructive. In addition, the dynamics of these quantities 
oscillate with time t when the initial state is |4 f i(0)) and the 
amplitudes of C, QD and GMD decay to a stable value after 
a long time evolution with the enhancement of 7, while that 
of CC, on the contrary, is enhanced with the increase of 7. 
And CC is maximal all the time but other quantities dissipate 
degressively and disappear eventually with t when the initial 
state is chosen as an alternative Bell state | "1/2 (0)) . Moreover, 
there is not a definite ordering of various quantities in thermal 
equilibrium, whereas there is a descending order of CC, C, 
GMD and QD under the intrinsic decoherence with a nonnull 
B when the initial state is chosen as 1^2(0)). 



V. CONCLUSION 

In summary, we have investigated various correlations mea- 
sured by C, CC, QD and GMD in a two-qubit Heisenberg 
XX Z spin chain in the presence of external magnetic field 
and DM anisotropic antisymmetric interaction both in ther- 
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